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We reconsider the original model of quadratic hybrid inflation in light of the WMAP three-year 
results and study the possibility of obtaining a spectral index of primordial density perturbations, 
n s , smaller than one from this model. The original hybrid inflation model naturally predicts n s > 1 
in the false vacuum dominated regime but it is also possible to have n 3 < 1 when the quadratic 
term dominates. We therefore investigate whether there is also an intermediate regime compatible 
with the latest constraints, where the scalar field value during the last 50 e-folds of inflation is less 
than the Planck scale. 
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I. INTRODUCTION 

The results from the WMAP three-year data [l[ have 
provided the first indication that the spectral index of 
primordial density perturbations, n Sf is smaller than one. 
The original hybrid inflation model [2j naturally predicts 
n s > 1 in the false vacuum dominated regime but it is 
also possible to have n s < 1 when the quadratic term 
dominates We are therefore interested in whether 
there is also an intermediate regime compatible with the 
latest constraints. Hybrid inflation models are attractive 
from a theoretical point of view because the inflaton field 
(j> is far below the Plank scale during inflation, in the false 
vacuum dominated regime, which makes these models 
easier to implement within supergravity Q . Therefore we 
also investigate whether there is an intermediate regime 
where the scalar field value during the last 50 e-folds of 
inflation is less than the Planck scale. 

In this paper we study the original hybrid inflation 
model with a quadratic potential for the inflaton. Wc 
want to analyze carefully this model in light of the new 
WMAP results, exploring the space of parameters of the 
model. We have run a code to calculate numerically the 
spectral tilt of the scalar curvature perturbations, n s , 
the running of the spectral tilt, a s , and the tensor-to- 
scalar ratio, r, assuming slow-roll, for a specific region of 
parameters of the model. We study the cases for which 
we have n s < 1 and compare the results obtained with 
the WMAP three-year results. We use the 68% and 95% 
confidence level contours from WMAP only and WMAP 
+ SDSS taken from 5]. Recently similar works have 
used the WMAP three-year results to put constraints on 
hybrid inflation models [1, 0, ■ 

We also study the inverted hybrid inflation model 
@ with a quadratic potential for the inflaton (see Ap- 
pendix) . 



II. THE HYBRID INFLATION MODEL 

The potential for the hybrid inflation model is given 

by 1 
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V(<t>, x) = -Ax 2 - Xlf + ^rn 2 ^ 2 + -A> 2 X 



(1) 



where (f> is the inflaton, \ is called the "waterfall" field, A 
and A' are coupling constants and xo and m are constant 
masses. 

It is assumed that x stays at the origin, which cor- 
responds to a false vacuum, while <fi rolls down from a 
initially large (positive) value until it reaches a critical 
value, 4> c — ■J A/ X'M, after which x becomes unsta- 
ble (the effective mass-squared becomes negative) and 
rapidly rolls down towards one of the true minima at 
X = ±Xo- Then <f> goes to zero and starts to oscillate 
while x will reach the true minimum. Inflation will end 
either with the instability or because of the end of slow- 
roll, as in the single field case, depending on which occurs 
first. 

Before <fr reaches <p c we can write the potential as a 
function of cf> only, 



V{<j>) = Af 4 + V 



(2) 



where we wrote the false vacuum energy density ^A%g as 
just Af 4 . We define the ratio between the false vacuum 
energy density and the inflaton energy density as 



E{<£) = - 



M 4 



m 2 (j) 2 



(3) 



So if E > 1 we have false vacuum dominated inflation 
and if E <C 1 we have almost chaotic inflation [lfj. 



A. The dynamics 



The dynamics of hybrid inflation are given by the equa- 
tion of motion of the inflaton (we are assuming that the 
"waterfall" field stays at the origin, so it does not evolve) 
and the Fricdmann equation, 
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-V'(<j>), 
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where H = a/a is the Hubble parameter, a is the scale 
factor, mpi is the Planck mass (which we set equal to 1), 
a dot represents a derivative with respect to time and a 
prime denotes a derivative with respect to the field <j>. We 
will use the slow- roll approximation for our calculations, 
which is given by the conditions 



the details of reheating at the end of inflation. Given 
that (j) c can be made arbitrarily small by suitable choice 
of A/ A' we leave </>5o, the value of <j> when cosmological 
scales leave the Hubble scale, as a free parameter to be 
determined by observations, subject only to the restric- 
tion that we must have at least 50 e-folds between (^50 
and (f>e- If 4><i does not exist then we can always have at 
least 50 e-folds and any value of a^o is allowed. 
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With these approximations, Eqs. (H} and © can be writ- 
ten as 



B. The perturbations 

The power spectrum of scalar curvature perturbations 
at horizon crossing (when the comoving scale k equals 
the Hubble radius, k = aH, during inflation), which is 
conserved on large scales in single field inflation, is g iven 
by, to leading order in the slow-roll parameters, [ll| 
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Then we can write, for the number of e-folds of expansion 
between two field values 4>i and 02, 
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For the potential in Eq. ([2]) we have 
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When the end of slow-roll occurs before the inflaton 
reaches <p c we identify the end of inflation with the con- 
dition e = 1, which occurs for the value of the inflaton 
field i 



(14) 



There is a second root at a smaller (f>, below which e 
becomes smaller than unity again, but numerically it has 
been found that slow-roll is not re-established before </> = 
Qi. 

We see that if 327rM 4 / 'm Pl m 2 > 1 then <j) t does not 
exist at all, so in this case inflation has to end by insta- 
bility. If 32n M 4 /m Pl m 2 ^ 1 it will end by instability 
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or by the end of slow-roll when 



when (j) = 

4> = 4>e if 4>t > 4> c - 

We assume that cosmological scales exit the Hubble 
scale at least 50 e-folds before the end of inflation. In 
practice the actual number of e-folds is dependent upon 



(15) 



where the subscript * indicates that the quantity is to 
be evaluated at horizon crossing. By virtue of the slow- 
roll conditions this formula gives a value of V-n which is 
nearly independent of k. For the potential in Eq. |2]), 
assuming slow-roll, Eq. (|15p can be written as, using 
Eqs. ©, ©, ©, © and (HH), 
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(16) 



So for each values of M and m we can find the values of 
</>50 (there are several possible values for each combina- 
tion of M and m) which satisfy the density perturbation 
amplitude, for scales of cosmological interest, from the 
WMAP three- year results [l[ (assuming that the running 
of the spectral tilt is zero) . 

The spectral tilt for the scalar curvature perturbations, 
the running of the spectral tilt and the tensor-to-scalar 
ratio can be written as, to leading order in the slow-roll 
parameters, [ll| 



n s = 1 - 6e* + 277*, 
a s = 16e*77* - 24e 2 + 2^ 2 , 
r = 16e„, 



(17) 
(18) 
(19) 



where £ is a higher order slow-roll parameter and is equal 
to zero for our specific potential, Eq. ([2]). Then for this 
potential we can write 
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C. Results and discussion 

We have run a code to calculate numerically 4>so j E and 
the parameters n s , r and a s at Hubble crossing (when 
<p = ), assuming slow- roll and using the definitions 
and results from the previous sections, for m between 
1(T 4 and 1(T 8 and M between 10~ 2 and 10~ 5 . For each 
values of m and M we have selected the value of (^50 such 
that the amplitude of the scalar curvature perturbations 
Eq. ([16]) obeys the WMAP three-year results We dis- 
card the cases for which the number of e-folds between 
4>5o and <f) £ is smaller than 50 and those with a blue spec- 
trum, i.e., we require n s < 1. The energy density ratio, 
Eq. (J3]), is evaluated at 5O . 

We first note that as one goes from larger to smaller 
values of M the density of points selected by the code 
gets larger because the selection of the values of M is 
logarithmic, as one can see looking to the plot in Figure 

m 



For the range of parameters for which [n s , r) is inside the 
95% confidence level contour we have 

0-95 < n s < 1 and < r < 0.3. (26) 
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FIG. 1: Contour plot of n a as a function of m and M for 0.9 < 
n 3 < 1, for the hybrid inflation model. The red, green, blue, 
yellow and light blue contours (from left to right) represent, 
respectively, n s equal to 0.91, 0.93, 0.95, 0.97 and 0.99. 

Analyzing the top plot in Fig. [5] we see that there is a 
small range of parameters for which (n s ,r) is inside the 
68% confidence level contour from WMAP only For 
this range we find 

0.9525 < n s < 0.9975 and 0.05 < r < 0.35. (23) 

For the range of parameters for which (n s , r) is inside the 
95% confidence level contour we find 

0-94 < n s < 1 and < r < 0.55. (24) 

Considering the results from WMAP + SDSS @ we see 
that, looking to the bottom plot in Fig. [21 the range of 
parameters for which (n s , r) is inside the 68% confidence 
level contour is smaller than in the previous case. For 
this range we find 

0-96 <n s < 0.995 and 0.05 < r < 0.175. (25) 



2 . 5 

2 

1.5 



1 




0.9 0.92 0.94 0.96 0.9S 



FIG. 2: Plots of r as a function of n s for 0.9 < n s < 1 (the 
black dots represent the different values of M and m), for the 
hybrid inflation model, with the 68% (yellow) and 95% (blue) 
confidence level contours from WMAP (top) and WMAP + 
SDSS (bottom), taken from 

Observing the plot in Figure [3] we see that the closer 
one gets to large values of M the larger are the values 
of r (i.e., the energy scale of inflation increases when M 
increases) and a s . Analyzing the plot in Figure [4] we 
see that the closer one gets to small values of </>5o the 
larger are the values of r (and so also the values of a s 
get larger). We also note that </>so is never much smaller 
than 1, i.e., the Planck mass. Looking to the plot in 
Figure [5] we note that the energy density ratio is never 
much larger than 1, which means that there is never a real 
false vacuum domination of the energy density, otherwise 
we would have n s > 1. We also see that the closer one 
gets to the false vacuum domination cases (which occur 
for large values of M) the larger are the values of r and 
the larger (and positive) are the values of the the running 

From this analysis we see that when we have false vac- 
uum domination and small values of 4>50t i- e -> when we 
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FIG. 3: Plot of r as a function of a s for 0.945 < n s < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the hybrid inflation model. 



FIG. 5: Plot of r as a function of E for 0.945 < n s < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the hybrid inflation model. 
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FIG. 4: Plot of r as a function of <f> 50 for 0.945 < n 3 < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the hybrid inflation model. 



are more distant from chaotic inflation, then the values 
of the the running a$o are large and positive, which is 
in contradiction with the WMAP three-year results ||. 
Also the values of the tensor-to-scalar ratio r tend to be 
very large for these cases, which is also in disagreement 
with the observations. 

When M goes to very small values we recover the re- 
sults for chaotic inflation, which correspond to the lower 
limit for r along the values of n s in Fig. O For these 
cases, for which 0so is very large and the energy den- 
sity ratio very small, we have small values for a^o and r, 
which is in very good agreement with the WMAP three- 
year results [H . 



III. CONCLUSIONS 

We have found that there is an intermediate regime 
for the original hybrid inflation model compatible with 
n s < 1, but as one approaches the false vacuum domi- 
nated limit within this regime the tensor-to-scalar ratio, 



r, and the running of the spectral tilt, a s , become large 
(with a positive running) , which is in contradiction with 
the WMAP three- year results This agrees with the 
results in [1, 0, § ■ We found a lower bound on the al- 
lowed values of r, which might be an observational signal 
for hybrid inflation because if there is an upper bound on 
the spectral index, n s < 1 , then we find a lower bound on 
the tensor-to-scalar ratio. This lower bound corresponds 
to chaotic inflation, which is in good agreement with the 
new WMAP results @. We also saw that it is difficult to 
get smaller than the Planck scale because decreasing 
050 requires us to increase the vacuum energy scale, M, 
and this increases the tensor-to-scalar ratio, r. 

The inverted hybrid inflation model (see Appendix) is 
in much better agreement with the WMAP three-year re- 
sults |5| than the original hybrid inflation model when we 
consider the false vacuum dominated regimes, especially 
because there is no lower bound on the allowed values 
of r. Moreover, in the inverted hybrid model there is no 
problem obtaining small values for c/j^q <ti 1 , which makes 
these model easier to implement within supergravity 0]. 
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APPENDIX: THE INVERTED HYBRID 
INFLATION MODEL 

If we invert the sign of the inflaton energy density in 
Eq. (J2J we get 



V((j>) = M 4 - im 2 2 , 



(27) 
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which corresponds to the potential for the inverted hy- 
brid inflation model with a quadratic potential for the 
inflaton. A particular case of this model are the small- 
field inflation models [ll[ , for which the potential in Eq. 
(|27p corresponds to a Taylor expansion about the origin 
and higher order terms are required to provide a poten- 
tial minimum at some cf> ^ 0, as required to connect to a 
reheating stage (which is necessary to make the model re- 
alistic). So the inflaton starts with an initial small (pos- 
itive) value and rolls down towards the minimum at a 
larger value, then it starts to oscillate and inflation ends 
with the end of slow-roll. 

All the equations obtained before for the original hy- 
brid inflation model, except Eq. (fl4| . are still valid but 
with to 2 replaced by —to 2 . To find the value of <fi e we 
search for the solutions of Eq. (fT2")) with e equal to f 
(there are two different solutions for this equation but 
the results do not depend on which solution we choose). 

1. Results and discussion 

In this case we have run the same code (with the ap- 
propriate changes in the equations) but for to between 
1(T 4 and 10~ 10 and M between 10~ 2 and 10~ 6 . The en- 
ergy density ratio, Eq. ([3]), is evaluated 50 e- folds after 
the inflaton has the value 

Note that also in this case the selection of the values 
of M is logarithmic, as one can see looking to the plot in 
Figure [5] 
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FIG. 6: Contour plot of n s as a function of m and M for 
0.9 < n s < 1 , for the inverted hybrid inflation model. The 
red, green, blue, yellow and light blue contours (from left to 
right) represent, respectively, n„ equal to 0.91, 0.93, 0.95, 0.97 
and 0.99. 

Analyzing the top plot in Fig. [7] we see that there is a 
much wider range of parameters for which (n s , r) is inside 
the 68% confidence level contour from WMAP only [5] 
than in the original hybrid case, especially because there 



is no lower bound on the allowed values of r, in contrast 
with the original hybrid inflation model. For this range 
we find 

0-935 < n s < 0.9875 and < r < 0.0775. (28) 

For the range of parameters for which (n s , r) is inside the 
95% confidence level contour we get 

0.9175 < n s < 1 and < r < 0.0775. (29) 

Observing the bottom plot in Fig. [7] we can see that 
there is also a wide range of parameters for which (n s , r) 
is inside the 68% confidence level contour from WMAP 
+ SDSS 0]. For this range we get 

0.935 < n s < 0.9925 and < r < 0.0775, (30) 

and for the range of parameters for which (n s , r) is inside 
the 95% confidence level contour we find 

0-92 < n s < 1 and < r < 0.0775. (31) 
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FIG. 7: Plots of r as a function of n s for 0.9 < n s < 1 (the 
black dots represent the different values of M and m), for the 
inverted hybrid inflation model, with the 68% (yellow) and 
95% (blue) confidence level contours from WMAP (top) and 
WMAP + SDSS (bottom), taken from [|. 
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Looking to the plot in Figure [8] we see that the closer 
one gets to large values of M the smaller are the values of 
r and a s (although we note that they are always small, 
even in the cases for which M is small). Analyzing the 
plot in Figure [9] we see that the closer one gets to small 
values of ^50 the smaller are the values of r (and so also 
the values of a s get smaller). We note that here ^50 
can be much smaller than the Planck mass. Observing 
the plot in Figure [TOl we see that the closer one gets to 
the more false vacuum domination cases (which occur for 
large values of M) the smaller are the values of r (and 
so also the values of a s ). In this case we note that the 
energy density ratio can be very large, so we can have a 
strong vacuum domination. 
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1, which means that there is a large false vacuum domi- 
nation of the energy density. Therefore we can conclude 
that this model is in much better agreement with the 
WMAP results than the original hybrid inflation model 
when we consider the false vacuum domination regimes, 
especially because there is no lower bound on the allowed 
values of r. Moreover, in this model there is no problem 
obtaining small values for </>5o <C 1, which makes these 
model easier to implement within supergravity 
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FIG. 10: Plot of r as a function of E for 0.945 < n s < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the inverted hybrid inflation model. 



FIG. 8: Plot of r as a function of a a for 0.945 < n a < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the inverted hybrid inflation model. 
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FIG. 9: Plot of r as a function of <j> 50 for 0.945 < n a < 0.955 
(red), 0.965 < n s < 0.975 (green) and 0.985 < n s < 0.995 
(blue), for the inverted hybrid inflation model. 

From the previous analysis we see that for the inverted 
hybrid inflation model we can have, specially for large 
M, very small values for r and a s , which is in very good 
agreement with the WMAP three-year results 0. For 
these cases ^50 is also very small (it can be much smaller 
than 1, i.e., the Planck mass) and E is much larger than 
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